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Abstract 



The analytical solution of the recently proposed ideal chain poly- 
mer mean-spherical approximation (Yu.Kalyuzhnyi, Mol.Phys., 94, 
735(1998)) is presented for the multicomponent mixture of charged hard- 
sphere linear chain flexible molecules. The solution apply to any mixture 
of chain molecules with arbitrary distribution of the charge and size of 
the beads along the molecular backbone. Closed form analytical ex- 
pressions for the internal energy, Helmholtz free energy, chemical poten- 
tials and pressure are derived. By way of illustration thermodynamical 
properties of several versions of the fluid of charged chain molecules of 
different length, including the molecules with uniform, diblock and al- 
ternating distribution of the charge are studied. Theoretical predictions 
are in reasonable agreement with available computer simulation predic- 
tions. We present also the liquid-gas phase diagrams for the systems 
with diblock and alternating distribution of the charge. 
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I. INTRODUCTION 



In recent years substantial amout of efforts has been focused on the development of the- 
oretical approaches predicting the equilibrium structure and thermodynamical properties of 
charged polymer solutions. This information enables one to esteblish the stability limits and 
to build the phase diagrams of such fundamentally and technologically importent systems as 
polyelectrolyte solutions (for example DNA, polyacrilic acid) and polyampholyte solutions 
(for example proteins). One of the simplest continuum models of charged polymer solutions 
is represented by the fluid of charged hard-sphere linear chain flexible molecules. In spite 
of its simplicity this model incorporates several essential features of charged polymer fluids, 
such as excluded volume effects, chain connectivity and flexibility, distribution of the charge 
along the chain backbone and Coulomb interaction. Recently several off-lattice theories for 
the fluid of flexible charged chain molecules, based on the integral-equation techniques de- 
veloped for the fluids of small molecules, have been proposed. These include extension |T[ of 
the polymer reference interaction site model (RISM) integral equation theory of Schweizer 
and Curro [Q], extension [Q,|J of the simple interpolation scheme (SIS) of Stell and Zhou || 
and extensions [p|-|T0|j of the multidensity integral equation theory for associating fluids of 
Wertheim [fUH . The latter studies are based on the analytical solution of the mean spher- 
ical approximation (MSA) version of the multidensity theory, the so-called polymer MSA 
(PMSA) |6|-p|,p!2| . In the case of uncharged system direct application of Wertheim's original 
formulation by Chang and Sandler gave rise to a polydisperse fluid of chains with 

a prescribe mean number of beads. A subsequent extension by Kalyuzhnyi and Cummings 
fT5|| yields a model that polymerizes to a fluid of chains of fixed length in the complete 
association limit. In that limit, two-density Wertheim's Ornstein-Zernike (OZ) equation be- 
comes identical to the Rossky-Chiles |]16[ version of the "proper" site-site integral equations 
fT7|| ("proper- RISM") as first noted by Kalyuzhnyi and Stell [|18[] and later discussed in more 
details by Stell and by Kalyuzhnyi and Cummings P0| . More recently complete associ- 
ation version of this extension was combined with the ideal chain approximation |TT| , [T3"||n 



3 



and generalized for the systems with arbitrary long-range potential outside the hard-core 
region |J . As a result the corresponding version of the PMSA was formulated and H0ye-Stell 
scheme [22] of calculating the simple fluids thermodynamical properties was extended in the 



frames of the PMSA M. Subsequently von Solms and Chiew |TTJ combine PMSA and ideal 



chain approximation in a way similar to that proposed by Chang and Sandler [£L3[]. The 
resulting polydisperse chain mixture was used to describe thermodynamical properties of 
the two-component mixture of fixed-length chain molecules with equally charged beads and 
oppositelly charged monomer counterions. Results of the theory appeare to be in reason- 



able agreement with computer simulation results [23|. We note in passing, that ideal chain 
approximation utilized by von Solms and Chiew []1D1 is somewhat different from the original 
Tl|JT3|| . In the corresponding Wertheim's OZ equation the authors take into account all 



one 



partial correlation functions, while due the original ideal chain approximation |TT| , [i3|j corre- 
lations involving at least one doubly bonded particle are neglected. More recently general 
solution of the PMSA |7]||] for polymerizing charged hard spheres, supplemented by the 
ideal chain approximation, was used to derive closed form analytical expressions for thermo- 
dynamical properties of charged chain fluid with a certain restrictions imposed on the sizes 
of the chain beads [p4 |. In the latter study the authors propose an improved version of the 



ideal chain approximation, which satisfies the Debye-Hiickel limiting law. 

In this paper we present an analytical solution of the ideal chain PMSA || for the mul- 
ticomponent mixture of charged hard-sphere flexible chain molecules. Unlike the previous 
studies proposed version of the PMSA is quite general and apply for any mixture of the 
chain molecules with arbitrary distribution of the charge and size of the beads along the 
molecular backbone. Our solution represents a complete association limit of the general 
solution obtained earlier |7],^] for the fluid of polymerizing charged hard spheres. Similar, 
as in the case of primitive electrolyte model |25[] this solution reduces to the solution of 
only one algebraic equation for the MSA-like screening parameter T. We also generalize re- 
cently developed method of calculating PMSA thermodynamics via the energy route |J for 
the multicomponent case and derived closed form analytical expressions for internal energy, 



Helmholtz free energy, chemical potentials and pressure in terms of the present solution. By 
way of illustration we consider several different versions of the fluid of charged hard-sphere 
chain molecules, including the molecules with uniform, diblock and alternating distribution 
of the charge along the chain backbone. 

The paper is orgenized as follows. In the next section we discuss detailes of the model 
to be studied and formulate ideal chain PMSA theory. Section III gives solution of the 
ideal chain PMSA and Section IV contains expressions for thermodynamic and structure 
properties of the system. Numerical results and discussion can be found in Section V and 
in Section VI we collect our conclusions. 



II. THE MODEL AND IDEAL CHAIN PMSA THEORY 

The model fluid we consider consists of the M-component mixture of freely jointed 
tangent hard-sphere chain molecules with each molecule of species a represented by m a 
charged hard-sphere sites of sizes a a a and charges ez a a . The total number density of the 
system is p = J2 a Pa, where p a is the number density of species a, and we assume charged 
neutrality conditions 

Ep^ = o (1) 

aa 

We denote the molecular species by the small letters a, b, c, ... taking the values 1,2,...,M 
and the site type in a given molecule a by a small greek letters a, (3, 7, which take the 
values 1,2, ...,m a . Thus the site has two indices, one denoting the species of the molecule 
and the other the site type in the molecule. For example the hard-sphere diameter of the 
site of type a belonging to the molecule of the species a is denoted by o a a . 

The site-site pair potential U^(r) between the sites of the type a and (3 belonging to 
the molecules of species a and b can be written in the form 

U$(r) = U^ )a \r) + U^ a \r) (2) 
where U^ ab (r) is the hard-sphere potential and ujfp ab (r) is the Coulomb potential 
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p 2 a b 

e Q r 



(3) 



with eo being the dielectric constant of the continuum. 

In this study we are using the multidensity version of the mean spherical approximation 
(MSA) [n]ri2|||fj (or polymer MSA (PMSA)) supplemented by the so-called ideal chain 
approximation P^JTBU^TI Pd . Both PMSA and the ideal chain approximation have been dis- 
cussed at length earlier and therefore we will omit any details and present here only the final 
expressions. The theory consists of the Ornstein-Zernike-like integral equation 



c 7 

and PMSA boundary conditions 
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where h^(k) and c^Jk) are the matrices with the elements being the Fourier transforms 
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of the elements of the matrices h^(r) and c^(r) 
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tij, ex. and E are the following matrices 
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Here the lower indices i,j each take the values 0, A, and B, and denote the bonding states 
of the corresponding particle, i.e. denotes the unbonded state and A and B denote the 
A-bonded and S-bonded states, respectively. The total partial correlation functions h^p. (r) 
are related to the site-site total correlation functions h^g(r) by the following relation 

O r ) = EW0 (6) 

ij 

The set of the OZ equations @ together with the PMSA closure conditions (|5]) represent 
our ideal chain PMSA theory for the fluid of charged linear chain hard-sphere molecules. 



III. SOLUTION OF THE IDEAL CHAIN PMSA FOR THE FLUID OF 
CHARGED CHAIN MOLECULES 



Recently a general solution of the PMSA for the multicomponent mixture of polymerizing 
charged hard spheres was published J7|||. Here we elaborate on this general solution by 
utilizing additionally the ideal chain approximation [TT| , |TB| J1 an d specializing it for the case 
of the model at hand. 

The general solution was obtained using Baxter's tecnique, which factorizes the initial 
OZ equation into a set of two equations 



Sj(|r|) = Qj(r) - EPcE / dr'QZJr')cx [Q 
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(M) = Q$(r) + 5>c£ / dr'J™(\r'-r\)cxQ%(r') 
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where T denotes the transpose matrix. The projections Srg(r) and J^( r ) 
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satisfy the following boundary conditions 

J^(r) = vrr 2 E + J 



Sj(r) 
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obtained from (|). Here = J^g(0) and the limit \i — > is to be taken at the end of the 
calculations. 

From the analysis of the equations (0) (|§D we get the following expression for the factor 
function Q^g(r) 



Q&{r) = 9 (o% - r) [q$(r) + t J] - (z*) T a 6 



(11) 



where 



/2 and the function q^(r) is defined in the range A^g < r < cr^ by 



= ^E T A^ (r - a%) (r - A&) + 0$ (r - <) . 



(12) 



Here E, and are the row vectors 



e = (1,0,0), K = K,< A ,< B ), s£ = «,o,o). 



Coefficients of the Baxter g-function can be expressed in terms of only one parameter, 
i.e. the screening PMSA parameter T |25|,[7]||. We have 
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Here T satisfies the following equation 
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where X^ is the row vector = ( X% , X aA , X aB ) . 

Thus solution of the PMSA for the present model reduces to the solution of the algebraic 



equation (21) for the screening parameter T. 
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IV. STRUCTURE AND THERMODYNAMICS 



A. Structure properties 



The contact values of the regular part of the partial pair distribution functions 
9a'i/3j(, 'af3+) follow from the relation (|8|) after differentiating it with respect to r and taking 

ab 
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The values of g^ if3 .(r) for r > c™^ can be calculated using the following relation 



(24) 



where ujfp ab (k) is the Fourier transform of the Coulomb potential U < fp ab (r) and expression 
for the function ^^.(k) = h^p (k) — c^p^k) follows from the set of equations (0) and (|J), 
written in the Fourier /c-space 
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Here p, h(k), c(k) and bQ(k) are the matrices with the elements p ab i p ] = ^ab^ap^ijPa-, (k), 
Ca^ik) and 
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respectively. Substituting expression for the fuction Q^p.( r ) (0) m t° ^ ne right hand side 
of the Fourier transformation (|27j) we find that 
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This result allow us to derive the expression for the structure factor 



Using the second of the relations ( |25| ) we have 
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B. Thermodynamic properties 



Expression for the excess internal energy AE follows from the standard relation 
AE 

= 27T 2^ PaPb 
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which after some calculations gives 
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where a a a N a a = £f =0 - < 
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The other thermodynamical properties of the system in question can be calculated fol- 
lowing the method developed earlier ||. This method generalizes H0ye-Stell MSA energy 
route to thermodynamics in the frames of the present PMSA approach. Originally this 
generalization was developed for the one-component case. Extension of the method in the 
multicomponent case is rather straightforward and yields the following expressions for the 
Helmholtz free energy A, pressure P and chemical potential \x a a in excess to their reference 
system values 
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where the quantities with the superscript (ref) denote the reference system quantities and 
the quantities with superscript (0) denote the corresponding zero charge PMSA quantities 
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The reference system is represented by the multicomponent mixture of uncharged hard- 
sphere chain molecules. Expressions (B3), Q3"5]) and (RoT) are quite general and apply for any 
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type of the potential outside the hard-core. The quantities, which enter these expressions 
can be written in terms of the present solution of PMSA for the model at hand. We have 
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Here expression ([41]) is obtained by differentiating of the equation (§) twice and considering 
the limit of r — > and expression (12) follows from the first of the equations fl2"5|) at 

k -> 0. 



V. RESULTS AND DISCUSSION 

To illustrate the above solution of the PMSA in this section we present the numerical 
results for thermodynamical properties of the several versions of charged chain fluid model. 
We consider two-component mixture of chain polyions with the beads of equal charge z v a = 
z p = — 1 and oppositely charged counterions z c a — z c — 1 (model Ml) and two versions 
of three-component mixture of chain polyions with diblock and alternating distribution of 
oppositelly charged beads, z v + = —z v _ = 1, and two types of counterions with opposite 
charges z\ = —zt = 1 (models M2 and M3). The densities of the counterions of these two 
models were choosed to be p\ = pi = \m p p v . Here the indices p and c denote polyion and 
counterion, respectively. We also consider the one-component version of charged chain fluid 
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models with diblock and alternating distribution of the charge along the chain backbone. In 
all cases the hard-sphere diameters of the chain beads a v a and counterions a c a were choosed 
to be equal, o v a = <j c a = a = 1. Schematic representation of the models studied are shown 
in figure 1. Thermodynamical properties of the uncharged versions of these models were 



calculated using TPTD theory of Chang and Sandler |14 



In figure 2 we compare osmotic coefficient results obtained using the present version 



of the PMSA, von Solms and Chiew version of PMSA [1(| and molecular dynamic (MD) 



computer simulation method |23| for the model Ml at a different chain length (m p = 16, 32 
and 64) and Bjerrum length Xb = (3* = 0.833. In the region of the intermediate densities 
our theory underestimate and the theory of von Solms and Chiew overestimate the value 
of the osmotic coefficient <fi = (3P/ p, where p = J2 a Pa- In the region of the low densities 
both theories predict for the osmotic coefficient its ideal gas value of = 0, while MD 
simulation gives the value, which is slightly above 1. This disagreement could be due to 
the uncertainties of the MD simulation in the diluted region. In general predictions of the 
present approach are more accurate than those of von Solms and Chiew approach Jl0| . 



In figures 3-5 we compare the density dependence of the osmotic pressure and osmotic 
coefficient at different values of Bjerrum length (Xb = 0.833,2.499) for different polyelec- 
trolyte models. In the case of low Bjerrum length (Xb = 0.833) one can see the linear 
dependence of the osmotic pressure with respect to the density on the log — log scale (figures 
3a-5a). In this region the osmotic pressure is almost independent of the chain length and 
distribution of the charge along the chain. With the increase of (Xb = 2.499) there is 
a departure from the linearity, which is substantial for the model Ml, slightly smaller for 
the model M2 and almost negligible for the model M3. At the same time one can see the 
substantial chain length dependence of the osmotic pressure in the case of the model M2. 
This dependence becomes smaller in the case of the model Ml and is absent for the model 
M3 (figures 3b-5b). Thus for the systems with more random distribution of the charge along 
the chain backbone the chain length dependence is smaller. With the increase of the chain 
length this dependence becomes weaker and it is reasonable to expect, that for m p > 64 it 
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will be negligible for all the models studied. 

In figur 6 the liquid-gas phase diagrams for the one-component version of the models 
with diblock and alternating distribution of the charge are presented. We consider the 
systems with molecular chain length m p = 2,4,6,8,10. Computer simulation results are 
available only in the case of the chains length m p = 2 |[26|| . In general theoretical results 
are in qualitative agreement with computer simulation results with certain disagreement for 
the slope steepness of the liquid branches of the coexistence curves; the slope of theoretical 
curve is steeper than that of the computer simulation curve. In addition, the position of the 
theoretical critical point is shifted towards lower values of the density and higher values of the 
temperature. According to the previous discussion the systems with diblock distribution of 
the charge have larger degree of nonideality in comparison with the systems with alternating 
distribution of the charge. This is reflected in the figure 6, where the critical temperature of 
the systems with diblock chains is higher than that of the systems with alternating chains 
of the same length. The corresponding values of the critical densities is larger in the case 
of systems with alternating chains. In both cases increase of the chain length leads to the 
increase of the critical temperature and slight decrease (increase) of the critical density for 
diblock (alternating) distribution of the charge. We note in passing, that due to the lack of 



the solution of the equation for screening parameter T (|21~D at lower densities we were not 
able to determine the coexistance curve for the systems with alternating chain charge of the 
chain length m p = 8, 10 below the temperatures shown in figure 6b. 

Finally in figure 7 we present the liquid-gas phase diagram for the system of trimers with 
the middle bead charge z\ — —2zf = —2z%, where z\ and z\ are the charges of the terminal 
beads, z\ = z\ = \. Similar as in the case of alternating charged chain system with m p = 8, 9 
there is no solution of the equation ( |2"T| ) in the low density region, therefore we were not 
able to calculate the coexistence curve for the temperatures lower than those shown in the 
figure. This model can be seen as a complete association limit of the primitive 2:1 electrolyte 



model. Recent computer simulation studies [^,|6|,[^|g9[] of the electrolyte primitive model 



suggest that in the vicinity of the coexistance region the system is highly associated and 
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fraction of charged clusters, including free ions, is negligible. Therefore the equilibrium 
properties of the system are determined by the presence of the neutral ionic clusters, to 
which the lower order clusters give the main contribution. Assuming that in the case of the 
electrolyte restricted primitive model (RPM) all such clusters are represented by the neutral 
ionic pairs the RPM liquid-gas phase diagram has been successfully reproduced by the 
Monte Carlo (MC) computer simulation method |26j and theoretically ||. Making similar 
assumption in the case of the 2:1 primitive electrolyte model the corresponding liquid-gas 
phase diagram can be modeled using the phase diagram obtained for trimers. Unfortunately 
computer simulation predictions for the entire phase diagram of the 2:1 electrolyte model 
currently is not available. Recently Camp and Patey [^9j present MC estimate for the critical 
temperature, which appeares to be T* = 0.1 ±0.01. The rough estimate of the ranges for the 
critical density can be obtained from the figure 5 of p9[ ; our estimate is p* = 0.115 ± 0.085. 
This estimate for the critical point together with Debye-Hiickel (DH) theory estimate |29| 
and MSA phase diagram are shown in figure 7. In addition, for the sake of completeness, 
we present the liquid-gas phase diagram for the RPM of electrolyte, obtained from the MC 
simulation method |27],[3(]], MSA and our PMSA for neutral diatomics |J. One can see 
that predictions of the PMSA is much better in comparison with usual MSA. In both cases 
DH predictions for the critical temperature is close to those of the PMSA, however the 
corresponding predictions for the critical density are much worse. 



VI. CONCLUDING REMARKS 

In this study we obtain an analytical solution of the ideal chain PMSA for the multi- 
component mixture of charged hard-sphere flexible linear chain molecules. This solution 
apply to any mixture of chain molecules with arbitrary distribution of the charge and size 
of the beads along the molecular backbone. We present closed form analytical expressions 
for thermodynamical and structure properties of the system. These expressions are used to 
calculate thermodynamics of several different versions of charged hard-sphere chain model. 
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We consider two- component mixture of chains with equally charged beads and oppositely 
charged counterions, two versions of the three-component mixture of chains with diblock 
and alternating distribution of the charge and two types of oppositely charged counterions. 
Comparison of the theoretical predictions for the osmotic pressure of the two-component 
model with corresponding computer simulation predictions shows that the present version 
of the ideal chain PMSA theory is reasonably accurate and gives somewhat better agreement 



than that of the theory of von Solms and Chiew [fTOj] . Effects due to the difference in Bjerrum 
length, molecular chain length and distribution of the charge was also studied. The largest 
deviation from ideal behaviour was demonstrated by the system with uniform distribution of 
the charge; behaviour of the system with alternating distribution of the charge shows much 
less nonideality. This conclusion is varyfied by the comparison of the liquid-gas phase dia- 
grams calculated for the one-component version of the models with diblock and alternating 
distribution of the charge; the latter model has substatially lower critical temperature than 
that of the former. With the increase of the chain length the critical temperature for both 
models increases, while the corresponding critical density slightly increases for alternating 
charge model and decreases for diblock charge model. Finally it was demonstrated, that 
liquid-gas phase diagram of the neutral trimers with doubly charged middle bead and singly 
charged terminal beads of the opposite sighn can be used to correct MSA results for the 
phase diagram of 2:1 primitive electrolyte model. 
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CAPTIONS TO THE FIGURES 



Figure 1. Schematic diagrams of the models with uniform distribution of the charge 
(Ml), diblock distribution of the charge (M2) and alternating distribution of the charge. 

Figure 2. Osmotic coefficient <fi = PV/NkT as a function of the packing fraction 
7] = 7rpcr 3 /6 for the model Ml at A# = 0.833 for m p = 16 (a), m p = 32 (b) and m p = 64 
(c). Predictions of the present theory (solid lines), von Solms and Chiew theory JTD| (dashed 
lines) and computer simulation predictions ]23| (diamons). 

Figure 3. Osmotic pressure P* = f3P a 3 (a) and osmotic coefficient = PV/NkT 
(b) as a function of packing fraction rj = 7ip T a 3 /6 (pT = J2a m aPa fc> r the model Ml from 
the present theory at A# = 2.499 (solid lines), \b = 0.833 (dashed lines) and \b = 
(dashed-dotted lines). For each set of lines from the top to the bottom at log 77 = —2 
m p = 8, 16, 32, 64. Symbols are computer simulation predictions ]23| for Xb = 0.833 and for 
m p = 16 (diamonds), m p = 32 (circles) and m p = 64 (squares). 

Figure 4. The same as in figure 3 for the model M2. 

Figure 5. The same as in figure 3 for the model M3. 

Figure 6. Liquid-gas phase diagram in T* = 1/(3* versus reduced density of the beads 
pi = m p pa 3 coordinates for one component chain fluids with diblock distribution of the 
charge (a) and with alternating distribution of the charge (b). From the top to the bottom 
at pi = 0.05 m p = 10,8,6,4,2. Solid lines are predictions from the present theory and 
diamonds are computer simulation predictions |26j for m p = 2. 

Figure 7 Liquid-gas phase diagram in T* = 1/(3* versus reduced ionic density p* = 
(p + +p_)cr 3 coordinates for restricted primitive model (RPM) of electrolyte (lower portions of 
the figure for T* > 0.08) and for the primitive model (PM) of electrolyte with 2:1 assymetry 
in charge (upper portions of the figure for T* < 0.08). Present theory (solid lines), MSA 
(dashed lines), Debye-Hiickel theory [^| (squares), computer simulation estimate for the 
critical point of 1:2 electrolyte PM |?9| (open diamonds) and for the phase diagram of 
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electrolyte RPM (open circles |27]], solid diamonds 
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FIGURES 




Figure 1 (Kalyuzhnyi and Cummings) 
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Figure 2 (Kalyuzhnyi and Cummings) 
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Figure 3 (Kalyuzhnyi and Cummings) 
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Figure 4 (Kalyuzhnyi and Cummings) 
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Figure 5 (Kalyuzhnyi and Cummings) 
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Figure 6 (Kalyuzhnyi and Cummings) 
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Figure 7 (Kalyuzhnyi and Cummings) 
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